home *** CD-ROM | disk | FTP | other *** search
/ Fritz: All Fritz / All Fritz.zip / All Fritz / FILES / PROGMISC / PCSSP.LZH / PC-SSP.ZIP / STATRAND.ZIP / GAUSS.FOR < prev    next >
Text File  |  1985-11-29  |  2KB  |  48 lines
  1. C
  2. C     ..................................................................
  3. C
  4. C        SUBROUTINE GAUSS
  5. C
  6. C        PURPOSE
  7. C           COMPUTES A NORMALLY DISTRIBUTED RANDOM NUMBER WITH A GIVEN
  8. C           MEAN AND STANDARD DEVIATION
  9. C
  10. C        USAGE
  11. C           CALL GAUSS(IX,S,AM,V)
  12. C
  13. C        DESCRIPTION OF PARAMETERS
  14. C           IX -IX MUST CONTAIN AN ODD INTEGER NUMBER WITH NINE OR
  15. C               LESS DIGITS ON THE FIRST ENTRY TO GAUSS. THEREAFTER
  16. C               IT WILL CONTAIN A UNIFORMLY DISTRIBUTED INTEGER RANDOM
  17. C               NUMBER GENERATED BY THE SUBROUTINE FOR USE ON THE NEXT
  18. C               ENTRY TO THE SUBROUTINE.
  19. C           S  -THE DESIRED STANDARD DEVIATION OF THE NORMAL
  20. C               DISTRIBUTION.
  21. C           AM -THE DESIRED MEAN OF THE NORMAL DISTRIBUTION
  22. C           V  -THE VALUE OF THE COMPUTED NORMAL RANDOM VARIABLE
  23. C
  24. C        REMARKS
  25. C           THIS SUBROUTINE USES RANDU WHICH IS MACHINE SPECIFIC
  26. C
  27. C        SUBROUTINES AND FUNCTION SUBPROGRAMS REQUIRED
  28. C           RANDU
  29. C
  30. C        METHOD
  31. C           USES 12 UNIFORM RANDOM NUMBERS TO COMPUTE NORMAL RANDOM
  32. C           NUMBERS BY CENTRAL LIMIT THEOREM. THE RESULT IS THEN
  33. C           ADJUSTED TO MATCH THE GIVEN MEAN AND STANDARD DEVIATION.
  34. C           THE UNIFORM RANDOM NUMBERS COMPUTED WITHIN THE SUBROUTINE
  35. C           ARE FOUND BY THE POWER RESIDUE METHOD.
  36. C
  37. C     ..................................................................
  38. C
  39.       SUBROUTINE GAUSS(IX,S,AM,V)
  40.       A=0.0
  41.       DO 50 I=1,12
  42.       CALL RANDU(IX,IY,Y)
  43.       IX=IY
  44.    50 A=A+Y
  45.       V=(A-6.0)*S+AM
  46.       RETURN
  47.       END
  48.